###########################################
# Doses exposure
###########################################

rm (list = ls())

library(Hmisc)
library(readstata13)
library(foreign)
library(lmtest)
library(sandwich)
library(stargazer)
library(msm)
library(ggplot2)
library(did)
library(lfe)
library(gridExtra)

# load data
load("one_week_protests_riots_2024March1.RData")

##############
# Table A10
##############

m1 = felm(support_eviction  ~ t_con_quintile_f + year_f + t_con_quintile_f*year_f + woman + age + education | 0 | 0 | idencuesta, data = d, keepCX = TRUE)
t1 = m1
m1 = tidy(m1,conf.int = T)
print(m1, n = 50)
